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THE GLOBAL RADIATION BALANCE OF THE EARTH 
ATMOSPHERE SYSTEM OBTAINED FROM RADIATION 
DATA OF THE METEOROLOGICAL 
SATELLITE NIMBUS II 

Ehrhard Raschke* and Musa Pasternak 
Laboratory for Atmospheric and Biological Sciences 


ABSTRACT 

The meteorological satellite Nimbus II measured the reflected 
solar radiation in the spectral range from 0.2 to 4. 0 microns and 
the emitted thermal radiation in the spectral range from 5. 0 to 30. 0 
microns over the entire globe during a period of about 2.5 months 
( 16 May - 28 July, 1966) . These measurements were obtained with 
a medium resolution scanning radiometer. Thus, to compute the 
outgoing fluxes of reflected solar radiation and of emitted thermal 
radiation, empirical integration methods had to be applied to the 
measured data. For the computation of the emitted thermal radia- 
tion a method previously developed by Wark et al. [23] was used. 
Functional relations were derived from balloon, airplane, and satel- 
lite measurements of reflected solar radiation to account for the 
angular dependence of the reflectivity of the earth-atmosphere sys- 
tem in calculating the flux of reflected solar radiation. Results 
were obtained for five subperiods each of half a month's duration. 

In this paper global maps are presented and discussed showing 
the geographic distribution of the albedo, the absorbed solar radia- 
tion, the emitted thermal radiation, and the radiation balance, ob- 
tained from the Nimbus II measurements during the subperiod 1 - 
15 June, 1966. All other results are discussed elsewhere [18]. 

The global average of the radiation balance for all five sub- 
periods was between -0.007 and +0.002 cal cm’ 2 min" 1 indicating 
that the earth- atmosphere system is nearly in radiative equilib- 
rium, These values are, however, within the limits of uncertainty 
of the computed fluxes of incoming and reflected solar radiation 
and the thermal radiation. The global average of the albedo was 
found to be about 30%, suggesting that other authors who found a 
value of about 35% from climatological data might have overesti- 
mated the cloudiness, especially in tropical latitudes. 


*On leave from the University of Munich as a Postdoctoral Resident Research 
Associate of the National Academy of Sciences. 
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THE GLOBAL RADIATION BALANCE OF THE EARTH 
ATMOSPHERE SYSTEM OBTAINED FROM RADIATION 
DATA OF THE METEOROLOGICAL 
SATELLITE NIMBUS II 


INTRODUCTION 

The earth-atmosphere system both receives energy from and looses energy 
into space by radiation. Solar electromagnetic radiation in the spectral range 
from 0. X -4.0 microns forms the major energy source. About 60 - 70% of it 
is absorbed in atmospheric layers below about 150 -200 km and at the ground. The 
rest is reflected and scattered back to space. Solar particle radiation influences 
highly the properties of atmospheric layers above 200 - 250 km. Its contribu- 
tion to the total radiation influx and also the contribution of radiation from other 
extraterrestrial sources, however, is small enough to be neglected in these 
considerations. Energy, is lost to space by emission of thermal radiation from 
ground and cloud surfaces and from several atmospheric gases. The main en- 
ergy transfer by thermal radiation takes place in the spectral range from 5 to 
50 microns. 

The absorption and reflection (and scattering) of incident solar radiation 
and the emission of thermal radiation are primarily effective at the ground and 
in the lower atmosphere up to about 40 km. Very shortwave solar radiation at 
wavelengths less than 0. 3 microns is absorbed in upper atmospheric layers, 
where it is responsible for many processes. It only amounts to about 1. 5% of 
the total incoming radiation flux. Thus the radiation balance Q at the top of the 
atmosphere, which is the difference between incoming and outgoing radiation, 
determines in the first order the radiation balance of the earth-atmosphere sys- 
tem below the 40 km level. It can be expressed by: 

Q(A.,0, d) = S - E(\,0,d) (1) 

In Eq. (1) S, R, and E are the fluxes of incoming and reflected solar radia- 
tion and of the emitted thermal radiation, respectively, crossing a horizontal 
element of area outside the atmosphere, whose geographic longitude and lati- 
tude are k and 4> respectively. The letter d designates the day of the year for 
which the balance computations are made, during which the sun has a certain 
declination and distance from the earth. 

Earlier computations of Q (e.g. see ref. [5], [lO], [l4]) were made with 
climatological data of various atmospheric parameters. All of these studies 
concluded that in general there is an excess of incoming radiation at low 
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latitudes and a deficit at high latitudes, the amplitude and location of which 
change with the seasons. This causes an energy transport by the general circu- 
lation in the atmosphere and by ocean currents. The yearly global average was 
found to be zero, suggesting that the earth-atmosphere system is in energy 
equilibrium over a period of a number of years. 

These computations, however, were based on observational material whose 
reliability especially over polar regions and over oceanic areas at that time was 
much in doubt. Therefore, satellite experiments were designed to measure the 
fluxes of reflected solar radiation and of emitted longwave radiation over the 
entire earth. Prior to the Nimbus II medium resolution radiation experiment, 
radiation experiments were flown on Tiros satellites which observed only parts 
of the earth’s surface between about 60° N and 60° S. Results obtained from 
these measurements (e.g. see ref. [ 3 ], [ll], [19]) only confirmed the magnitude 
of the earlier computed values of the outgoing longwave radiation. The outgoing 
flux of reflected and scattered solar radiation was computed from the satellite 
data assuming isotropy of the reflection properties of surfaces and the atmos- 
phere above them (Lambert surface) . The results so obtained were too low to 
balance on a global scale the incoming solar radiation (S) with the outgoing ra- 
diation fluxes (R and E) . Since there was no calibration on board the satellite, 
it was suspected that the instruments suffered a degradation in their response. 
Based on this and on the assumption of an annual global radiation balance Q = 0, 
Bandeen et al. [s] derived correction factors with which the satellite measure- 
ments of outgoing solar radiation were multiplied to determine the total outgoing 
flux of solar radiation from satellite measurements. Rasool and Prabhakara 
[l9] and Arking [l] have treated Tiros radiation data in a similar way. 

Theoretical computations (e.g. see ref. [7], [9]) as well as some measure- 
ments from balloons or airplanes ([4], [6], [2l]) however, indicate a predomi- 
nant dependence of the reflection properties of surfaces and of the atmosphere 
above on the zenith angle and azimuth angle (with respect to the sun) of obser- 
vation, and also on the elevation of the sun above the horizon. Arking [l] con- 
firmed these findings with a statistical analysis of Tiros IV measurements of 
the reflected solar radiation in the spectral interval from 0.2 - 4. 0 microns. 

Thus, it is suggested, that the assumption of a Lambert-surface for calcu- 
lations of the outgoing reflected flux of solar radiation is the reason for the ap- 
parently low results. As a consequence, for the computations of the outgoing 
flux of solar radiation from Nimbus n measurements efforts have been made to 
take into account the angular dependence of the reflection properties of the 
earth- atmosphere system. 

The medium resolution radiometer in Nimbus II provided measurements of 
the emitted longwave radiation and reflected solar radiation over the entire 
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globe during a period of about 2. 5 months (16 May - 28 July, 1966), This period 
covers the beginning and development of the summer and winter seasons in the 
Northern and Southern hemispheres, respectively. But it is too short, to draw 
any conclusions about the annual radiation balance. 

This entire period was divided into five subperiods, each of half a month’s 
duration, for which global distributions of S, E, E, and Q (Eq. 1) and some re- 
lated quantities were computed. In this paper only results which were obtained 
from the Nimbus II measurements of the period 1 June - 15 June, 1966 are pre- 
sented and discussed. All other results are presented in more detail in a com- 
prehensive report on this project [18], 


METHODS OF COMPUTATION 

Nimbus II was launched on May 15, 1966 into a sunsynchronous , nearly 
polar and circular orbit. At a mean height of 1140 km above the surface of the 
earth, it crossed the equator in the northward direction at around local noon 
and in the southward direction at around local midnight. Due to its orbital peri- 
od of 107 minutes and due to its height, generally each area element on the earth 
was observed at least once in the daytime and once in the nighttime within a 24 
hour interval, finis within 24 hours reflected solar radiation was observed at 
least once over each area of the illuminated portion of the earth. Emitted long- 
wave radiation could be measured from the same area both during day and night. 
The Nimbus II experiment is described in more detail elsewhere [17]. 

Two channels of the five-channel medium resolution radiometer were de- 
signed to measure the reflected solar radiation (0.2 - 4. 0 microns) and the 
emitted longwave radiation (5. 0 - 30. 0 microns). The instantaneous angular 
field of view of the radiometer subtends an overall angle of about 4 degrees. 
Therefore, the computation of both fluxes through an area outside the atmos- 
phere, which is perpendicular to the earth’s radius vector through it, would 
require measurements from ail directions to perform the integration over the 
entire hemisphere, as expressed mathematically by: 

r 2 tt r n/2 

Flux = / / N(0,t/O sin 6 cos 6 66 6 ^ (2) 


In Eq. 2 6 is the zenith angle of observation at the observed area, \p is the azi- 
muth counted from a fixed position (e.g. the horizontal projection of the solar 
ray incident upon the observed point) , and N(0,t/>) stands for the total, radiance 
of reflected sunlight or of emitted thermal radiation determined from the satel- 
lite measurements. 
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An integration according to Eq. (2) at the height of the satellite would cause 
a great loss of detail in all measured radiation fields (e.g. , see [15]). There- 
fore, it was assumed here, that the reference level is located on the surface of 
a spherical earth with a radius of 6371 km. This includes the assumption that 
this surface has the same physical properties (reflection, scattering, refraction, 
emission) as the earth-atmosphere system. Now the outgoing fluxes of reflected 
solar radiation (R) and of emitted longwave radiation (E) were considered to be 
radiation fluxes from an element of area into the upward hemisphere above it, 
The geographic coordinates of this element of area are determined from a know- 
ledge of the motion of the spacecraft in orbit and of the phase of the scanner 
mirror in the reduction of the original analog signals. The computation methods 
are only very briefly outlined here; they are discussed in more detail elsewhere 
[18]. 


The outgoing flux of emitted longwave radiation has been computed by the 
method of Work et al* [23], This method consists of two major steps: 

1, From the radiance measured in the spectral range from 5. 0 to 30. 0 
micron, the radiance over the entire spectrum is computed. For this a relation 
between both quantities has been derived from computations of both the "filtered" 
and "unfiltered" radiance emerging from a set of more than 100 atmospheric 
profiles into different directions {0 ) . Some results of these computations , which 
were carried out by Lienesch [l3], are presented as points {0 ~ 0) and as en- 
circled points {d~ 78. 5°) in Fig. 1. These points approach narrowly a straight 
line because the range of sensitivity of the filter (5,0 - 30. 0 microns) covers 
that spectral region within which about 90% of the totally emitted longwave ra- 
diation emerges to space. 

2, The second step includes the integration of radiances emitted into all 
directions. Here only the dependence on the zenith angle 6 was taken into ac- 
count by a function obtained statistically from Tiros VII measurements [13], 

No dependence of N was assumed with respect to the azimuth angle. T<; obtain 

a representative average of longwave radiation emitted during a 24 hour irtorval, 
the results from daytime and nighttime measurements were averaged, weighting 
each with the length of day and of night, respectively. 

The outgoing flux, of reflected and scattered solar radiation also was com- 
puted as an average for a 24 hour interval. The method used here consists of 
3 steps: 


1. Due to its narrow field of view, the radiometer in Nimbus II measured 
only that portion p of solar radiation incident on an observed area which is re- 
flected into the direction of the satellite. It is related to the measured radiance 
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# 


N f (#,■//) and to the solar radiation S f (d) incident in the same filter range (£) by 




Nf(^W0 


cos 


S' 


S f (d) 


(sr" 1 ) 


(3) 


Here £' is the zenith angle of the sun at time of measurement, and d expresses 
the day of measurement as described in the introduction above. S f has been 
obtained by integrating over all wavelengths the product of the spectral response 
of the instrument with Johnson's [l2] data on the extraterrestrial spectral solar 
irradiance. Since the spectral response of the instrument covers that range, 
within which about 99 % of the energy of solar radiation enters the atmosphere, 
it is assumed in all further steps, that p is a mean value for the entire solar 
spectrum. 


2. p is only the reflectance into one direction. The total reflectance r {p) 
is the ratio between the total reflected and scattered solar radiation in the up- 
ward hemisphere and the incident solar radiation. The functional relationship 
r (p) has been determined using conversion factors which were derived from 
many measurements of the angular dependence of the reflected solar radiation 
made from airplanes and balloons ([4j, [6], [2l]) and from Arking's [ 2 ] statisti- 
cal analyses of Tiros IV measurements. (For a Lambert surface r = 

3. In general, the r of an area element over a given underlying surface de- 
pends strongly on the zenith angle of the sun £. However, the functional rela- 
tionship between r and £ also depends upon the nature of the underlying surface, 
the turbidity of the air, etc. In these calculations only one function describing 
the dependence of r on the sun's zenith angle has been used. It was derived 
from the measurements and analyses mentioned above. Figure 2 compares our 
function with results by Arking and Levine, showing a general increase of total 
reflectance of the earth- atmosphez e system as the sun approaches the horizon. 
In this third step the total reflectance r was computed for all other zenith angles 
of the sun which are possible between sunrise and sunset over a certain area. 

The Luxes of incoming and of outgoing solar radiation were computed 
on the basis of the currently accepted value of the solar constant: 

S Q = 2. 0 cal cm" 2 min" 1 . Corrections were applied for the actual earth-sun 
distance on a given day d and for the effects of refraction, based upon a stand- 
ard atmosphere. 

The use of the above mentioned generalized functional relations undoubtedly 
causes errors where these functions do not coincide with the actual conditions . 
But, it is assumed that these errors are kept within an acceptable limit averag- 
ing all single measurements of a longer time interval within larger areas. 
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Therefore, a standardized grid program was used in which all single measure- 
ments are averaged within fields of the size of about 5 degrees of longitude and 
2 to 5 degrees of latitude. 


DISCUSSION OF THE RESULTS 


Albedo an d Absorbed Solar Radiation 

The albedo A of an area is defined here as that part of the daily incoming 
flux of solar radiation which is reflected back to space: 


A(\, c/>, d) 


d) 


100 [%] 


(4) 


It is characteristic for the reflection properties of the earth-atmosphere system 
at this area and determines the absorbed (ABS) part of the incoming flux by: 


ABS(\,<M) = ^1 - 

= S(\,4>,d) - R(^,d) 


(5a) 

(5b) 


Global maps of these quantities are shown in Figs. 3 and 4. 

It is known from other investigations that clouds, ice and snow surfaces are 
very good reflectors of incoming solar radiation. Therefore, the high albedoes 
in Fig. 3 can be interpreted to be caused either by clouds and/or by snow and 
ice surfaces. Cloud fields associated with the Intertropical Convergence Zone, 
the Monsoon over Southeast Asia, and the polar frontal zones of both hemis- 
pheres cause albedoes of more than 40%. Over the ice and snow covered polar 
cap north of 70° N, 60-80% of the incident radiation is reflected back to space. 
The reflectance of the surfaces here is increased by the low elevation of the sun 
throughout the entire day (Fig. 2) . The desert areas of North Africa, Arabia 
and also of Australia exhibit albedoes of more than 30-40%. 


Water surfaces have albedoes of less than 10%. Therefore, over the less 
cloud-covered subtropical oceans only 10-20% of the incident solar radiation is 
lost by reflection. The higher values for the albedoes of latitudes north of 50° N 
and south of 40° S indicate increased cloudiness. The general pattern of the 
albedo field in Fig. 3 and consequently of the field of absorbed solar radiation 
(Fig. 4) shows many irregularities over the Northern Hemisphere, due to the 
different in'” rences of sea and land surfaces on the circulation and therefore on 
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the reflection properties of the earth-atmosphere system. The preponderance 
of water in the Southern Hemisphere causes a more zonal pattern which is even 
more pronounced in the field of absorbed solar radiation (Fig, 4) , due to the 
strong meridional gradient of incident solar radiation at these latitudes. Since 
measurements obtained at solar zenith angles £* of more than 80° were omitted, 
both fields shew no data south of 58° S. The pola\ 5 cap south of 68° S was not il- 
luminated by the sun at this time of the year. 

Table 1 

24-IIour Average of the Incoming Flux of Solar Radiation (S) 

During the Period 1 - 15 June 1966 at Given Latitudes. 

Values in cal cm' 2 min* 1 . 



North 

South 

Latitude: 

o 

O 

80° 

O 

c 

O 

o 

20° 

0° 

20° 

O 

o 

60° 

68° 

Flux: 

0.75 

0.74 

0.69 

0.71 

0.67 

0.57 

0.41 

0.22 

0.06 

0 


If the albedo of the Northern Hemisphere were uniform, the North Pole would 
daily absorb during this time of the year even more solar energy than the sub- 
tropics (Table 1). The geographic distribution of the albedoes, however, causes 
a maximum of solar radiation of more than 0. 54 cal cm- 2 min- 1 to be absorbed 
over a latitude belt between about 10° and 35° N (Fig. 4). This belt is inter- 
rupted over North Africa and Arabia, and Southeast Asia where deserts and 
clouds, respectively, cause a high albedo. From this belt the amount of ab- 
sorbed solar radiation decreases toward both poles more or less continuously, 
due primarily to the increase of albedo (northward) and the decrease of the in- 
coming flux (southward) . 

In general meridional profiles of the absorbed solar radiation show that 
over the oceans more radiation is absorbed than over land at the same latitude. 

Outgoing Longwave Radiation 

The absorption spectrum of the earth’s atmosphere has a wide window be- 
tween 8 and 13 microns , through which radiation emitted from ground or cloud 
surfaces penetrates to space only weakly influenced by absorption and emission 
in intervening atmospheric layers. In this part of the spectrum also Planck's 
function of the spectral radiance of thermal radiation has its maximum for 
temperatures between 230° K and 300° K normally occurring in the atmosphere 
and at the ground. Therefore, the outgoing flux of longwave radiation can be 
considered to be in first order a function of the temperature of the underlying 
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surface. This surface might be a cloud, land, or ocean. This interpretation 
seems to be confirmed by a high correlation between simultaneous measure- 
ments of the upward going radiation in the spectral intervals 8-13 microns and 
7.5 - 32.5 microns, found by Moller and Has chke [l6] from radiation data of the 
satellite Tiros XXL Linear correlation coefficients of between 0. 85 and 0, 95 
were found for data of different orbits. Thus the pattern of the outgoing long- 
wave radiation as shown in Fig. 5 can be interpreted in terms of the tempera- 
ture of underlying cloud, land, or ocean surfaces. Pronounced mimima of out- 
going longwave radiation have been found over both polar caps, especially over 
the Plateau of Queen Maud Land. Very high and cold cloud fields associated 
with the intertropieal convergence zone and with the monsoon over Southeast 
Asia occur in Fig. 4 as mimima at low latitudes. The emitted flux is high over 
the sparsely cloud covered areas of the subtropical oceans and over the warm 
land surfaces at the same latitudes. 

Radiation Balance Q 

The net radiation flux at the top of the atmosphere or the radiation balance 
of the earth-atmosphere system can be found by subtracting the outgoing long- 
wave radiation of each area from the corresponding value of absorbed solar 
radiation. Figure 6 shows the geographic distribution found for the period 1 - 
15 June, 1966. 

These results show that all over the globe between 10° S and about 70° N the 
absorbed solar radiation exceeds the emitted longwave radiation. This leads to 
a surplus of energy which is partly stored in the oceans, ground, and atmos- 
phere, and partly transported by the general circulation and by ocean currents 
into regions of radiation deficit. Maxima of radiation excess are found over the 
subtropical oceans, while the hot deserts in the same latitudes show even a 
slight radiation deficit. This result is in good agreement with the work of 
Budyko [s] which showed an annual radiation deficit over the same regions. 

The highest radiation deficit of the Northern Hemisphere was found over 
Greeland, where only about 20-25% of incoming solar radiation is absorbed. 

At southern latitudes between 55° and 70° deficits of more than 0. 24 cal 
cm- 2 min -1 result from relatively high surface temperatures of ice and probab- 
ly open water and from the very low amount of incoming solar radiation. Since 
at this time of year the South Pole, especially the Plateau of Queen Maud Land, 
has very low surface temperatures of about 205° K, it shows a lower radiation 
deficit than the ice belt north of it. 

Meridional profiles of zonal averages of Q obtained from these computations 
and from the results of other authors are compared in Fig. 7. They show, that 
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the other authors predicted in their results very similar meridional profiles, 
but found lower values for the radiation excess between 0° and 30° N. This is 
probably caused by higher values of the albedo found by these authors over these 
latitudes (e. g. , London 31 - 34% vs Nimbus II 22 - 27%) . The global average of 
the radiation balance (Table 2) was found for all five subperiods to be between 
+0. 002 and -0. 007 cal cm -2 min* 1 , confirming that the earth atmosphere sys- 
tem is close to a radiation balance. These values, however, are differences 
between two very large numbers, the globally absorbed solar radiation and 
emitted longwave radiation, and are within the uncertainties of those large num- 
bers. Thus, no final conclusion on a net heating or cooling of the entire sys- 
tem during each of these subperiods can be drawn from these values. 

Table 2 

Global Radiation Balance from Nimbus II Measurements 



Albedo 

% 

Absorbed Solar 
Radiation 
cal cm* 2 min* 1 

Outgoing Longwave 
Radiation 
cal cm* 2 min* 1 

Radiation 
Balance 
cal cm* 2 min* 1 

May 16-31 

30.1 

0.341 

0.339 

+0 . 002 

June 1-15 

30.6 

0.337 

0.342 

-0.005 

June 16-30 

30.1 

0.338 

0.345 , 

-0.007 

July 1-15 

29.1 

0.343 

0.346 

-0.003 

July 16-28 

29.5 

0.342 

✓ 

0.345 

-0.003 


In Table 2 the values obtained for the global albedo are also listed. Ranging 
from 29. 1 to 30.6%, they are all lower than earlier estimates of the annual 
planetary albedo, which were found to be around 35% (e.g. , see ref. [8]). Direct 
comparisons of results obtained from the Nimbus II measurements with earlier 
estimates obtained from climatological data are only possible for the Northern 
Hemisphere, for which London [14] computed the radiation balance and some 
related quantities, for the summer season (June-August). This comparison is 
given in Table 3. 

The albedo value found by London is higher than those obtained from the 
Nimbus II measurements. Conversely, the total outgoing flux of longwave ra- 
diation as determined by London is lower than those obtained from the Nimbus II 
measurements. 

Clouds influence both the albedo and the outgoing longwave radiation. From 
the consideration made above, it might be suspected that surface observations 


9 



used by London may have overestimated the low latitude cloudiness, thus pro- 
ducing too high values of the albedo and too low values of the outgoing flux of 
longwave radiation. 


Table 3 

Radiation Balance of the Earth- Atmosphere 
System Over the Northern Hemisphere 




Albedo 

% 

Absorbed 
Solar 
Radiation 
cal cm* 2 min* 1 

Outgoing 
Longwave 
Radiation 
cal cm* 2 min* 1 

Radiation 
Balance 
cal cm* 2 min -1 

London 

June-July- 

August 

36.3 

0.411 

0.339 

+0.072 

Nimbus II 

May 16-31 

32.8 

0.443 

0.345 

+0.098 


June 1-15 

32.8 

0.452 

0.349 

+0.103 


June 16-30 

32.2 

0.458 

0.353 

+0.105 


July 1-15 

31.5 

0.458 

0.355 

+0.103 


July 16-28 

30.9 

0.446 

0.355 

+0.091 
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ALBEDO |%] 

NIMBUS 1 1 - 15 JUNE 1966 



Figure 3. The albedo of the earth-atmosphere system during 
the period 1-15 June, 1966. Values in %. 





ABSORBED SOLAR RADIATION [LY MIN" 1 ] 
NIMBUS 1 1 - 15 JUNE 1968 



Figure 4. Solar radiation absorbed daily in the earth-atmosphere system durinq the 

1 - 15 June, 1966. 










ET RADIATION FLUX AT THE 
OP OF THE ATMOSPHERE 



LATITUDE 

Figure 7. Zonal averages of the radiation balance of the earth-atmosphere system. 
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